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ABSTRACT 

From observations made with the ESPaDOnS spectropolarimeter, recently installed 
on the 3.6-m Canada-France-Hawaii Telescope, we report the discovery of a strong 
magnetic field in the Of?p spectrum variable HD 191612 - only the second known 
magnetic O star (following 9 1 Ori C). The stability of the observed Zeeman signature 
over four nights of observation, together with the non-rotational shape of line profiles, 
argue that the rotation period of HD 191612 is significantly longer than the 9-d value 
previously proposed. We suggest that the recently identified 538-d spectral-variability 
period is the rotation period, in which case the observed line-of-sight magnetic field 
of —220 ± 38 G implies a large-scale field (assumed dipolar) with a polar strength of 
about —1.5 kG. If confirmed, this scenario suggests that HD 191612 is, essentially, an 
evolved version of the near-ZAMS magnetic O star 8 1 Ori C, but with an even stronger 
field (about 15 kG at an age similar to that of 1 Ori C). We suggest that the rotation 
rate of HD 191612, which is exceptionally slow by accepted O-star standards, could 
be due to angular-momentum dissipation through a magnetically confined wind. 
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1 INTRODUCTION 

Magnetic fields of O stars ca n strongly imp act on the physics 
of th e stellar interiors (e.g., |S pruit 2003) an d atmospheres 
(e.g., iBabel fc Montmerld [1997^7 and hence on the stars' 
long-term evolution (e.g.. iMaeder fc Mevnetl l200l liool) . 
However, quantifying such effects critically depends on our 
knowledge of the properties of such fields. 

From a theoretical point of view, the origin of mag- 
netic fields in O stars is still rather uncertain. One possibil- 
ity is that they are fossil remnants of the star-formation 
stage, as propose d for magnet ic chemically-peculiar Ap 
and Bp stars (e.g. Mestel 119991) : if this were the case, we 
would expect to find about 10% of O stars with dipolar- 
type magnetic topologies with super-equipartition strengths 
of at least several hundred G. Another option, gaining 
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recent support among theoreticians, is that such fields 
are produ ced by dynamo processes, within the convec- 
tive cores JCharb onncau fc MacGrego3l200llh or m a pu- 
tative subsurface shear laye r jM^jcDroiaklfcMullanl 12004 
iMullan fc MacDonaldl l2005lh magnetic topologies would 
then depend on rotation rate and feature a significant 
toroidal component. 

From an observational point of view, however, very few 
direct constraints exist on the strength and topology of mag- 
netic fields in O stars, even though such fields are often 
invoked as a potential explanation for many otherwise enig- 
matic pheno mena (e.g., th e unex pected properties of the X- 
ray spectra, ICohen et alj (120031) '). While several attempts 
have been made to characterize the fields that could be 
hosted at the surfaces of bright, archetypal O stars such 
as £ Pup and ( Ori, no detections have yet been obtained 
with only one exce ption: the very young object 8 1 Ori C 
JPonati et al.l Eo02h One reason for this is that absorption 
lines of O stars are both relatively few in number in the op- 
tical, and generally rather broad (due to rotation or to some 
other type of as yet unknown macroscopic mechanism; e.g., 
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Table 1. Journal of observations. Columns 1-5 list the date, he- 
liocentric Julian date, UT time, exposure time, and peak signal to 
noise ratio (per 2.6 kms — 1 velocity bin) for each observation. The 
last column lists the rms noise level (relative to the unpolarized 
continuum level and per 41.6 kms" 1 velocity bin) in the circu- 
lar polarization profile produced by Least-Squares Deconvolution 
(Section . 



Table 2. Lines used for Least-Squares Deconvolution. The line 
depths (column 3) were directly measured from our spectra while 
the Lande factors (column 4) were derived assuming LS coupling. 
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iHowarth et"aTlll997l) . decreasing dramatically the size of the 
Zeeman signatures that their putative fields can induc e. 

The recent discovery bv lWalborn fc et all l|2004h that 
the Of?p star HD 191612 varies between spectral types 
06 and 08 with an apparently strict periodicity is sig- 
nificant in this context. The Ha and He I lines show a 
particular ly strong modulation , reminiscent of that seen in 
9 1 Ori C JStahl fc et al.]|l996l) ; HD 191612 thus appears at 
first glance as a very good candidate for the detection and 
investigation of a hot-star magnetic field. However, the de- 
rived period of the spectral variability, 538 d, is so long (by 
O-star standard s) that rotational modulation was consid- 
ered unlikely bv lWalborn fc et al.l Instead, they preferred 
to associate the spectroscopic variability with an eccentric 
binary orbit, hypothesizing that tidally induced non-radial 
pulsations (NRP) at periastron may induce periodically en- 
hanced mass loss. 

However, the lack of any clear orbital radial-velocity 
varia tions, down to a level of a few kms -1 flWalborn et alJ 
2003), is a significant hurdle to this model, especially since 
rather large eccentricity is implied by the NRP hypothesis. 
Thus, even though very slow rotation is at odds with conven- 
tional wisdom, the possibility that the observed variability 
may be related to rotational modulation of a magnetic star 
is by no means excluded. Motivated by its similarities to 
8 1 Ori C, we therefore included HD 191612 in a list of can- 
didates to search for magnetic fields in O-type stars using 
ESPaDOnS (Donati et al., in preparation), the new, high- 
efficiency spectropolarimeter recently installed on the 3.6-m 
Canada-France-Hawaii Telescope (CFHT). 

In this paper, we first present the observations and the 
Zeeman detection we obtained (Sec. [2j . We then review the 
constraints on rotation (Sec. perform simple modelling 
of our observations (Sec.0J, and discuss the implications of 
our results for rotational evolution of hot stars with strong 
radiation-driven winds (Sec.[SJ. 
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2 OBSERVATIONS 

Spectropolarimetric observations of HD 191612 were col- 
lected in 2005 June, as part of a four-night run aimed at 
investigating the magnetic fields of hot stars. (Results ob- 
tained on the other stars observed in the same run will be 
published separately.) At the time, the instrument was suf- 
fering a 1.3-mag light loss compared to the optimal perfor- 
mance obtained during the engineering runs (Donati et al., 
in prep.). This problem was not evident until the time of the 
run; now fixed, it turned out to be due to severe damage to 
the external jacket of optical fibres linking the polarimeter 
with the spectrograph (probably caused during movement 
of the instrument up and down from the Cassegrain focus). 
Current ESPaDOnS performance is therefore significantly 
better than the results given here imply. 

HD 191612 was observed each night; altogether, 13 
circular-polarization sequences, each consisting of 4 indi- 
vidual subexposures taken in different polarimeter config- 
urations, were obtained (see Donati et al., in prep., for 
details). All frame s were processed using Libre ESpRIT 
iDonati et al .119971: Donati et al., in prep.), a fully automatic 
reduction package installed at CFHT for optimal extraction 
of ESPaDOnS spectra. The peak signal-to-noise ratios per 
2.6 kms -1 velocity bin range from 200 to 450, depending 
mostly on weather conditions (see Table [T|. 

Least-Squares Deconvolution (LSD: lDonati et alJll997f) 
was applied to all observations. This requires the input of 
a line list, which was constructed manually to include the 
few moderate to strong absorption lines that are present in 
the spectrum of HD 191612. As, essentially, we aim to probe 
the photosphere of HD 191612, lines appearing in emission 
(such as the N ill lines at 463.41 and 464.06 nm) were not 
included, even though in some cases they may result from 
selective emission processes in or near the photosphere. The 
strong Balmer lines, all showing clear emission from the wind 
and/or circumstellar environment at the time of our observa- 
tions, were also excluded from the list, as were other features 
showing P-Cygni profiles (such as He n 468.57 nm and He I 
587.56 nm). This left only 12 usable spectral lines, whose 
characteristics are summarized in Tabled 

From those lines we produced a mean circular polariza- 
tion profile (LSD Stokes V profile) as well as a mean un- 
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Figure 1. LSD circularly-polarized and unpolarized profiles of 
HD 191612 (top, bottom curves respectively) on 2005 June 22- 
25. The mean polarization profile is expanded by a factor of 100 
and shifted upwards by 1.05 for display purposes. 



polarized profile (LSD Stokes I profile) for each spectrum. 
Given the breadth of the spectral lines of HD 191612 (av- 
erage full width at half depth is ~150 kms -1 ), the LSD 
profiles were produced on a spectral grid with a velocity bin 
of 41.6 kms -1 ; as this bin size is comparable to the typical 
broadening of intrinsic profiles from hot stars, no informa- 
tion from small-scale magnetic features potentially present 
at the surfac e of the star is expe cted to be lost in this binning 
process (e.g. iDonati et ai]Eo 02) . The resulting noise level in 
the LSD Stokes V profiles is about 0.02% of the unpolarized 
continuum per 41.6 kms -1 velocity bin (see Table0. 

Averaging together all LSD Stokes V profiles recorded 
on each of the four nights of observation (with weights pro- 
portional to the inverse variance of each profile) yields rela- 
tive noise levels of 2.3, 0.9, 1.1 and 1.2 respectively, in units 
of 10 _4 7 c . Compared to a null-field, V = profile, the re- 
sulting Stokes V signatures yield reduced-x 2 values of 0.9, 
1.9, 1.6 and 1.5, respectively, implying clear Zeeman detec- 
tions on each of the last three nights, though not the first. 
However, all four V profiles are mutually consistent. Averag- 
ing all V signatures yields a final noise level of 0.6 x 10 _4 7 c , 
and a reduced-x 2 of 2.7. We can thus safely conclude that 
the star is magnetic; the resulting average Stokes V profile 
is shown in Fig. and corresponds to a longitudinal field of 
-220 ± 38 G. 

Note that, even though we excluded lines showing an 
obvious contribution from the wind, the LSD Stokes I pro- 
file is still significantly asymmetric (with the blue wing ex- 
tending to higher velocities than the red one). Evidently, 
most 'photospheric' lines are likely to be affected to a small 
extent by the stellar wind. The Stokes V profile, showing a 
longer tail at blue wavelengths, is also slightly impacted by 
the same effect; however, adding or removing a few lines in 
the list does not change the overall shape of the detected 
Zeeman signatures (at the noise level), indicating that this 
perturbation has no significant consequences for the analysis 
presented here. 



3 THE ROTATION PERIOD OF HD 191612 

IWalborn fc et ah] i2004h inferred a rotation period for 
HD 191612 of ~9 d, from an estimate of the radius and 
a measurement of 77 kms -1 for the line- width parameter 
wsin(i). However, assuming that the magnetic topology of 
HD 191612 is essentially dipolar (as found for all other mag- 
netic OB stars identified to date), a rotation period of 9 d 
would imply significant evolution of the observed (projected) 
field configuration over a 4-d run, provided that the mag- 
netic dipole is not aligned with the rotation axis. We there- 
fore conclude that the rotation period of the star is most 
probably significantly longer than 9 d, consistent with the 
fact that variability is detected in neither our unpolarized 
spectra, nor in previous spectr oscopic data sets col lected 
over intervals of a few days (e.g. IWalborn et al.ll2003l) . 

This conclusion is further supported by the obvi- 
ously non-rotational shape of the LSD Stokes I profile of 
HD 191612 (Fig. 0; if 'turbulence', of some as yet unde- 
termined physical nature, dominates the line broadening, 
then clearly the line- width parameter usin(i) significantly 
overestimates t he true project ed equatorial rotation veloc- 
ity, n e sin(i) (cf. lHowarth 2004), leading to an underestimate 
of the rotation period. To examine this possibility further 
we compare the observed profile of C IV 580.1 nm to some 
simple models (Fig. |5J . This line is rather symmetrical, and 
is expected to form relatively deep in the atmosphere, so 
that it is probably more representative of the hydrostatic 
regions than i s the LSD I profile. W e took an OStar2002 
model profile llLanz fc Hubenvll2003l) . with modest scaling 
to match the observed line depth, and applied a rotational 
convolution with v c sin(i) = 77 kms -1 . The resulting pro- 
file does indeed match the observed line width quite well, 
but the overal match of the line shape is very poor. How- 
ever, a simple model of gaussian isotropic turbulence, with 
zero rotational broadening, while lacking any strong physi- 
cal justification, nonetheless provides an excellent match to 
the observations. 

Taken together, these arguments strongly suggest a ro- 
tation period significantly longer than the 9-d value pre- 
viously proposed. We are therefore tempted to explore the 
consequences of identifying the rotation period of HD 191612 
with the 538-d period recently identified for this star 
JWalb orn fc et al 1 l2004h . on which several spectroscopic in- 
dexes (and in particular Ha) are observed to vary. Although 
indicators of the implied very slow rotation are not yet 
overwhelming (e.g., constancy of the observed field could 
result from a dipole aligned with the rotation axis), the 
phenomenological similarities with the young O-type star 
6 1 Ori C, now strengthened by the magnetic field we have 
detected in HD 191612, indicate that the possibility of rota- 
tional modulation merits further consideration. Thus, while 
the proposal clearly needs further scrutiny (e.g., through the 
monitoring and potential detection of the rotational modu- 
lation of the Zeeman signature), we nevertheless now explore 
the consequences of a 538-d rotation period for HD 191612. 

4 MODELLING THE MAGNETIC FIELD OF 
HD 191612 

We propose that HD 191612, like 8 1 Ori C, hosts a magnetic 
field significantly tilted with respect to the rotation axis, and 
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Figure 2. The C IV 580. lnm profile, compared to models broad- 
ened by rotation (v c sin(i) = 77 kms -1 , dashed line) and by 
isotropic gaussian turbulence (ov = 45 kms -1 , solid line). The 
feature at —230 kms -1 is a diffuse interstellar band. 



that most of the observed variability results from the interac- 
tion of the radiative ly driven wind with the stellar magnetic 
field ( as described inlBabel fc Montmerleil997l:lDonati et alJ 
| 2002l:lud-Doula fc Owockill2002l : iTownsend fc Owocyi2005t 
iGaene et all2005l) . In this framework, the wind coming from 
each stellar hemisphere is deflected by the field towards the 
magnetic equator, where it produces a strong shock, an X- 
ray emitting post-shock region (reaching temperatures of 
10 6 to 10 7 K) , and a cooler and denser disk in the magnetic 
equator where the wind plasma piles up before being ejected 
away from (or accreted back onto) the star. In a generic way, 
we can then ascribe the observed Ha variations, in particu- 
lar, to the varying aspect of the cool circumstellar disk; the 
emission attains a maximum strength when the disk is seen 
pole-on, and a minimum when the disk is seen edge on. This 
could arise through recombination in a moderately optically 
thick disk. 

With this assumption, the longitudinal field of 
HD 191612 should be maximum at maximum Ha emis- 
sion, i.e., at phase 0.5 in the ephemeris obtained by 
IWalborn fc et al.J (120041) : JD min = 2, 448, 315 + 538S. More- 
over, the strongly reduced Ha emission observed around 
phase 0.0 indicates that the magnetic equator hosting the 
Ha emitting disk is seen close to edge-on in this viewing 
configuration. We therefore suggest that, as for 6 Ori C, 
the angle of the magnetic axis to the rotation axis /3 is close 
to 90° — i, where i is the angle of the rotation axis to the line 
of sight, to ensure that the magnetic equator is periodically 
seen edge on (at phase 0.0) by the observer. Given the large 
amplitude of the Ha variability, we can also conclude that 
neither i nor /3 is likely to be small. 

For an initial, schematic, modelling attempt, we there- 
fore propose for HD 191612 the simplest possible magnetic 
geometry, with i = f3 = 45° (as for 9 1 Ori C). In this con- 
text, the magnetic equator is seen edge on at phase 0.0, and 
the magnetic pole is f acing the observer at phase 0.5, in 
the lWalborn fc et al.J ll2004l) ephemeris. Our magnetic ob- 
servations were taken at phase 0.725; by fitting the detected 
Zeeman signature to a magnetic dipole model, whose single 
remaining free parameter is the field strength at the visible 
pole, we infer that the intensity of the dipole is — 1.5±0.2 kG. 

An alternative is to imagine that the rotation period 



is actually twice 538 d, so that the magnetic poles rotate 
alternately into the line of sight. However, we consider this 
configuration to be unlikely. First, it would require (3 to be 
very close to 90° , so that the disk is viewed in the same orien- 
tation (and thus produces the same amount of Ha emission) 
when each pole comes closest to the observer. Secondly, this 
option would produce edge-on disk viewing episodes that are 
much shorter than those obtained in the /3 — 90° — i case 
(where the disk approaches only asymptotically the edge-on 
configuration), in poor agreement with observations. 



5 DISCUSSION 

Several observational peculiarities of HD 191612 find a nat- 
ural ex planation in th e fram ework of this model. For in- 
stance, IWalborn et al.l l)2003t) estimate that the mass loss 
of HD 191612 is about 3 times stronger at phase 0.5 (i.e., 
when the magnetic pole and associated open field lines face 
the observer) than at phase 0.0 (i.e., when the magnetic 
equator and associated closed field lines are seen edge on). 
This is precisely what is expected in the cont ext of a mag- 
netica lly confined wind (see, e.g., Fig. 6 of iDonati et alJ 
( 2002 )). Similarly, the sh ape of the Hipparcos light curve 
( Walb orn fc et ah] i2004l) ) can be qualitatively explained by 
electron scattering in the disk, redirecting stellar photons to 
the observer at phase 0.5 and away from the line of sight at 
phase 0.0. 

Since a rotation period of 538 d is rather long by O- 
star standards, it naturally raises the question of whether 
the magnetic field is responsible for angular-momentum 
loss that produced the slow rotation. Focussing again on 
the very young star 8 1 Ori C suggests clues to answer 
this question. Since both stars have similar masses (of 
about 40 Mq), HD 191612 can, to first order, be con- 
sidered as an evolved version of 8 1 Ori C (whose age 
does not exceed 0.2 Myr). In this evolu tion, the radius in- 
creases from about 8 Ra (for 9 1 Ori C, iHowarth fc Pririial 
1 19891: IDonati et alJl2002l) to about 18 R Q (for HD 191612, 
Walb orr^^ll l^^^) while the temperature decreases from 
about 45,000 K iHowarth fc Priniall98sl:lDonati et all2002l) 
to about 35,000 K iWalborn et alJ 120031) : this is in rough 
agreement with evolutionary mo dels of massive stars (e.g. 
ISchaller et allll992l IClareti (20041) . from which we then de- 
rive an age of about 3-4 Myr for HD 191612. This is in good 
agreement with age estimates for the Cyg OB3 association 
( Masse^^t^dJ^y)95|, e.g.), of which HD 191612 is a member 
]Humphrevslll97Sl) . 

This scenario would imply that HD 191612 hosted a 
field of about 15 kG when on the main sequence at an age 
similar to that of 8 1 Ori C today. However, the correspond- 
ing change in the moment of inertia (about a factor of 3, 
taking into account the simultaneous evo lution of the frac- 
tional gyration radius k from 0.29 to 0.17;|ciarct 2004) does 
not of itself explain the change in rotation period between 
9 1 Ori C (15 d) and HD 191612 (538 d, if our model is 
correct) by more than an order of magnitude. 

To investigate whether the magnetic field may be re- 
sponsible for the required angular-momentum loss, we can 
evaluate the magnetic braking timescale through the simple 
expression 
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where _Ra is the Alfven radius (i.e., the distance up to which 
the wind is magnetically confined) and m is the effective 
mass-loss rate (determined by taking into account only the 
wind plasma that effectively leaves the star and thus con- 
tributes to angular- momentum loss). To evaluate the effect 
of the magnetic field on the wind of HD 191612, it is useful 
to consi der the wind magnetic confin ement parameter r/, de- 
fined bv lud-Doula fc Owo cki (2002) and characterizing the 
ratio between magnetic-field energy density and the kinetic- 
energy density of the wind: 

T)= B^RS/MVao, (2) 

where B oq is the equatorial magnetic field, M is the average 
mass loss rate, and v x is the terminal wind velocity. If r\ is 
significantly larger than 1, the magnetic field confines the 
wind within the Alfven rad ius, which roughly scales as r\ L ^ 
Jud-Doula & OwockilEool. 

While both HD 191612 and 6 1 Ori C have similar ter- 
minal wind velocities (vao ~ 2,500 kms -1 ), HD 191612 
exhibits significan tly stronger mass loss (by about an or - 
der of magnitude: iDonati et aT]l200l Iwalborn et al.ll2003t) . 
HD 191612 is also twice as large as, and features a ~ 1.5x 
stronger magnetic field than, 6 1 Ori C, implying that rj is 
similar for both stars (rj ~ 10), and th us that their Alfven 
radii are roughly e qual (at about 2 : IDonati et al.ll2002L 
iGaene et al J 120031 . It also implies that the effective mass- 
loss rate of HD 191612, m, corresponding to plasma evacu- 
ated through field lines opened by the wind (i.e. with a mag- 
netic colatitude smaller than about 45°), is about 25% of the 
actual surface mass flux, M, or ~ 1.5 x 10" 6 M yr" 1 . Given 
the age of HD 191612, the magnetic braking timescale we 
derive, of order 1 Myr, indicates that the field of HD 191612 
can, potentially, generate a strong enough brake to have 
slowed the star down to the currently observed rotation rate. 

One may argue that the angular-momentum loss of 
HD 191612 may not have been as strong as now throughout 
its past life. While still on the main sequence with a ra- 
dius and mass-loss rate comparable to those of 9 1 Ori C, 
HD 191612 would have had a magnetic field ~15 times 
larger than that of 1 Ori C, implying that its wind mag- 
netic confinement parameter was about r\ ~ 5000, and its 
Alfven radius reached about 15 ii* at that time. However, 
the correspondingly low effective mass-loss rate (of order 
10 -8 Mq yr -1 , or ~5% of the total surface mass flux) would 
still have imposed rotational braking on a timescale of order 
1 Myr, the larger Alfven radius roughly compensating for 
the smaller mass-loss rate. We therefore conclude that the 
angular-momentum loss of HD 191612 did not drastically 
change throughout the life of the star, and that HD 191612 
has had ample time to spin down since it was born. 

This interpretation does not explain why 6 1 Ori C itself 
is apparently rotating more slowly than normal O stars; it 
is far too young for its magnetic wind to have influenced its 
rotation rate significantly. We speculate that, at some stage 
in its formation process, the magnetic interaction between 
the forming star and its accretion disk may have prevented 
the star from accumulating as much angular momentum as 
normal, weakly magnetic, hot stars, in a mechanism s imilar 
to that proposed for magnetic Ap stars iStepienl2000l) . This 



speculation needs to be elaborated properly, however, and 
tested with adequate observations to see if it can realistically 
explain the slow rotation of newly born magnetic O stars. 

Further spectropolarimetric observations of HD 191612, 
sampling the 538-d period of spectrum variability, are ob- 
viously needed to establish firmly whether this is indeed 
the rotation period of HD 191612; to put quantitative con- 
straints on the magnetic-field geometry of this newly dis- 
covered magnetic hot star; and to test the preliminary con- 
clusions proposed in this paper. X-ray observations (such as 
those already undertaken by Naze et al., in prep.) over the 
538-d period will also help to constrain the magnetospheric 
physics and geometry, through the fluxes and spectral-line 
shapes formed in the postshock hot-plasma torus, as well 
as from the period ic eclipses of the torus that the star may 
naturally provide l|Babel fe Montmerlelll997t IDonati et alJ 
12001 iGagne et al]l2005T)~ 
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